Two-Dimensional t-J Model at Low Electron Density 
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The phase diagram of the two-dimensional t-J model is determined accurately at low electron 
density by a combination of analytic and numerical techniques. The ground state exhibits three 
phases in the limit of zero electron density: an unpaired state at small J/t, a gas of s-wave pairs 
for 2 < J/t < 3.4367, and a phase separated state at larger interaction strengths. Bound states of 
larger clusters are never realized, and the instabilities present at small densities are discussed. 



PACS numbers; 71.10.+X, 71.45. Gm, 74.20.-z 

Models of strongly-correlated electrons have received 
renewed attention since the discovery of the cuprate su- 
perconductors Q . However a reliable method of calcu- 
lating ground-state properties of such models on dense 
large-size systems has not yet been demonstrated. 

In this paper, we concentrate on the low electron den- 
sity limit of the two-dimensional t-J model, one of the 
simplest models with strong correlation ^-^. This limit 
appears far from the small hole densities thought to be 
applicable to cuprate superconductivity, but it is impor- 
tant to calculate the entire phase diagram of the model. 
A complete description of the low electron density limit, 
where much can be calculated analytically, yields a solid 
basis from which to move into the potentially more inter- 
esting regions of the model. Additionally, we introduce 
and test a numerical technique that we expect will prove 
to be a powerful tool for investigating strongly correlated 
models at higher electron densities. 
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acts only on the constrained Hilbert space with no dou- 
bly occupied sites. The operator c\^ creates an elec- 
tron of spin a on site z, S^ is the spin operator, and 
iii = X^cr '^L'^icr- ^^^ sum ovcr {ij) enumerates nearest- 
neighbor bonds of a square lattice. 

As J/t -^ 0, the electrons separate to avoid interac- 
tions in the limit of zero density, and the ground-state 
energy per particle is the kinetic energy of free electrons, 
£n=i = — 4i. Assuming a macroscopic number of elec- 
trons, the system phase separates completely in the large 
interaction strength limit, and the energy per particle, 
£n=oo = — 1.16934J± 0.00003J, comes purely from the 
interaction term in (|l|). This value is the result of the 
most accurate calculation of the ground-state energy of 
the Heisenbcrg model ^ shifted hy —jj per bond. 

For intermediate interaction strengths, the situation 
becomes more complicated. It is well known that for 
J/t > 2, pairs of electrons bind into an s-wave state 
P,H. However, previous work only cursorily examined 
the transition from this pairing phase to the completely 
phase-separated state at large J/t. 



This paper presents a thorough calculation of the phase 
diagram of the two-dimensional t-J model at zero elec- 
tron density. We calculate the energy per particle of all 
potential intermediate phases as a function of J/t. We 
show that no intervening phase is stabilized between the 
gas of s-wave pairs and the completely phase-separated 
state, and we determine the phase boundaries accurately. 
Finally we discuss possible instabilities at low densities. 

The energy of a bound electron pair may be found 
analytically, and we give the exact solutions below. In 
the limit of large J/t, a pair binds into a singlet dimer 
with an energy per particle of £n=2 — —\J ■ Since the 
chemical potential of a bound triplet is the same in the 
strong interaction limit, we dismiss three-particle bound 
states as a viable intermediate phase. 

The next competitive intervening phase to consider is a 
gas of bound quartets. Four electrons bind into a square 
configuration with an energy £n=4 — ~iJ P^r particle 
in the large J/t limit. Since this value is significantly 
lower than the pair chemical potential, quartets warrant 
careful scrutiny as a potential ground-state phase. We 
use a Green's Function Monte Carlo method to solve the 
four-body problem on finite lattices and extrapolate to 
infinite system size. Close to the critical J/t for four 
electrons to bind, the weakly bound quartets are very 
extended objects. An accurate calculation of this critical 
interaction strength on the infinite lattice requires pre- 
cise ground-state energies of four electrons on very large 
lattices. With the computational method described be- 
low, we can calculate ground-state properties on lattices 
with dimensions as large as 20 x 20. The results enable us 
to determine that bound quartets are never energetically 
stabilized in a macroscopic system. 

The paper is organized as follows: The solutions of two 
and four particles on the infinite lattice are outlined. We 
then calculate and discuss the low-density phase diagram. 

The equations of motion of two electrons yield two 
bound states for sufficiently large interaction strengths, 
one with s-wave symmetry and one with cia;2_j,2 symme- 
try 0,0-^. The analysis is simplified by defining the re- 
duced s-wave pair binding energy 6s = {—8t — Es)/8 > 0, 
where Es is the total energy of the s-wave pair. The re- 
duced d-wave pair binding energy Sd is defined similarly. 



On the infinite lattice, the s-wave solution satisfies 



1 
J 



1 



fAi^ 



(2) 



where K{k) is the complete elliptic integral of the first 
kind y. From now on, we express all energy scales in 
units of t. This expression may be expanded close to the 
critical J for pairs to bind. To order Sg, we find 
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where J* = 2 is the critical interaction strength to bind 
s-wave pairs. An expansion may also be made for large 
binding energy [||j|. 

After a similar analysis to that given in Ref. Q, the 
large-lattice limit of the d-wave solution reduces to 



1 
1 



2 + Sd^f VTTS, 



1 + 6d/2 



l + Sd 



(4) 



where E{k) is the complete elliptic integral of the second 
kind. An expansion for small binding energy yields 
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where the critical interaction strength for d-wave pairing 
is j;;' = 27r/(4-7r)«7.32. 

We solve the four-electron problem with a version 
of Green's Function Monte Carlo (GFMC) §-|l|l. We 
project a trial state onto the ground state by operating on 
it repeatedly with the Hamiltonian. This generates a se- 
ries of increasingly accurate approximants to the ground 
state labeled by integers \p) — HP\'^). 

It is useful to expand the trial state in terms of the 
exact eigenstates: 
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where j^I^o) is the ground state, |^i) is the first excited 
state, etc., and the a^'s are the expansion coefficients. 
Rewriting the projected states in this way, we see 
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where Ei is the energy of the i'th eigenstate. Thus \p) ap- 
proaches the lowest eigenstate for large p provided ao 7^ 
and |i?i>o| < 1^0 1 for all excited state energies ii^i>o, a 
condition satisfied by the t-J model for J > 0. 

Care is needed to choose a trial state having maximal 
overlap with the true ground state. We try to write a 
very general form for the trial state but restrict it to be 
a spin singlet with zero total momentum. We use the 
Jastrow-pairing state: 

i*> =n/(^«-^j-)^4nK+«k4^ciki)io) (8) 

i<j k 



FIG. 1. Plot of the GFMC output in the form of equation 
(td) for four electrons on a 20 x 20 lattice with J/t — 5.5. 
The data converge to the total ground-state energy at large 
power p. The solid line is the fit to ( |ll[ ) taking five terms 
in (h). The inset shows the scaling of the total ground-state 
energy with linear system size L. The fit to (113) is used to 
extrapolate to the infinite system. 



where c^^ is the Fermion creation operator and P4 
projects the state onto the subspace with four particles 
||l3| . We use mixed boundary conditions since the GFMC 
converges faster with a closed shell. 

The Jastrow factor /(r) correlates all pairs of parti- 
cles, yielding a total-spin-singlet state. We satisfy the t-J 
model's hard-core constraint by requiring /(r = 0) =0. 

The GFMC needs a positive-definite guiding function, 
which we take to be 



^' = Y[r{r.^r,)Y[s^ir.,^-~r,,^). 
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In general /^(r) 7^ /(r)- Since it is not important to 
guide with a spin-singlet function, we use the additional 
spin-dependent Jastrow factor s^ (r) in the guiding func- 
tion. The primed indices enumerate electrons of one spin. 

Auspicious choices of trial and guiding states can dras- 
tically reduce the statistical errors in the GFMC. We use 
much more general and potentially far superior states 
than previous work. Since the electrons on a given size 
lattice only occupy a finite number of locations r, we 
let the factors in dsh and (0) vary independently at each 
distance or wave vector not related by symmetry. We 
apply all rotation and mirror symmetries to the Jastrow 
factors, but only the mirror symmetries about the axes 
and parity to the Fermion pairing fields Uk and t;k- For 
example, the Fermion pairing fields may be any real lin- 
ear combination of an s and (1^2 _y2 pairing state, but 
time-reversal breaking and d^y symmetries are excluded. 

On a 20 X 20 lattice, each Jastrow factor has 400 pa- 
rameters, as do the Fermion pairing fields. The symmetry 
restrictions reduce the 800 total parameters in each state 
to 172 independent parameters in the trial state and 128 
in the guiding function. To optimize the parameters, we 
minimize the variance of the local energy |14| . For suf- 
ficiently large J/t, we find the spin-independent Jastrow 
factors in (H) and (^) bind the electrons in the optimized 
state, while the pairing fields Uk and Vk and the guiding 
function's spin-dependent Jastrow factor s^(r) provide 
the internal correlation of the bound quartet. 

The statistical errors increase exponentially with in- 
creasing power p, so we use the GFMC output at small 
powers to extrapolate to infinite power. We assume that 
the expansion m) can be approximated by a small num- 
ber of terms, typically about five, and we fit the output 
to determine the coefficients a^ and energies Ei. We in- 
clude enough terms in (0) to make the systematic error 
due to omitting terms smaller than the statistical errors. 



FIG. 2. Energy per particle for each phase at zero elec- 
tron density. The horizontal dotted line is the kinetic en- 
ergy of free electrons, and the dot-dashed line is the potential 
energy of the macroscopically phase-separated system. The 
short-dashed and long-dashed lines are the energies of s-wave 
and d-wave two-particle bound states, respectively. Their in- 
tersections with the free electron energy are marked with di- 
amonds, and the intersection of the s-wave pair energy with 
the phase-separated energy is marked with the square. The 
circles with error bars are the four-particle energies, and the 
solid line is the best fit of the quartet binding energies to ([14l) . 
The fit intersects the s-wave energy at the triangle. 



Previous work Mhof fit the GFMC output to 
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Fig. nl shows a nicely converged example in this form. 
While this function has the advantage of converging to 
the ground-state energy for large p, it is quite difficult to 
fit ( [ll|) with more than two terms in (^. We instead fit 
the simpler function 
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to determine all the coefficients a^ and energies Ei. 

For each value of J/t, we calculate the ground-state 
energy on square lattices of dimensions L x L, where L 
ranges from L = A to L = 20. We expect the energy 
of a bound state to converge exponentially with linear 
system size. To determine the infinite-lattice ground- 
state energy i5(L = oo), we fit the finite results to 



E{L) = E{L = oo) - Dcxp{-L/£,) 
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where E{L) is the ground-state energy on an ixi lattice, 
and D > and ^ > are parameters. The inset in Fig. 
|l| shows a sample extrapolation with system size. 

The energies of each phase are shown in Fig. g. The 
two-particle bound-state energies plotted are the exact 
expressions obtained by solving the integral equations (Q) 
and (0). The s-wave pair binds with a weak logarithm 
beyond the critical interaction strength J^ = 2, while 
the d-wave binding energy turns on with a healthy linear 
term for J > J^ Ri 7.32. The four-particle energies are 
extrapolated to infinite system size. Four-particle ener- 
gies for J < 5.21 are statistically indistinguishable from 
the s-wave pair energies and are omitted. 

Only three of the phases considered are realized as 
ground states. For J < 2, the electrons are free, minimiz- 
ing their kinetic energy. At larger interaction strengths, 
they bind into s-wave pairs, and finally the system under- 
goes macroscopic phase separation to form a Heisenberg 



cluster. The transition to complete phase separation oc- 
curs at JP* = 3.4367 ± 0.0001, which is determined by 
equating the phase separated energy with the solution 
of (H). The error of the transition point is due to the 
uncertainty in the best estimate of the ground state en- 
ergy of the Heisenberg model H. This value disagrees 
significantly with previously published results MM . 

It is obvious from Fig. that four-particle bound states 
are never stabilized in a macroscopic system. To deter- 
mine the critical interaction J^ for two pairs of electrons 
to bind into a four particle state, we fit the quartet bind- 
ing energies close to the transition to the form 
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where A > 0, -B, and C are parameters. The binding 
energy is given by (54 = —{Epf=4 — 2Es) > 0, where i?Af=4 
and Es are the total infinite-lattice ground-state energies 
of the four-particle and s-wave pair states, respectively. 
The form of the fitting function is a generalization of 
the first three terms in the expansion of the s-wave pair 
binding energy (ph, and includes the first two terms of 
the d-wave pair binding energy (H). 

The transition for pairs to bind into quartets takes 
place at J| = 5.07±0.15. The statistical error is large due 
to the singular nature of the fitting function (p4|). Still 
it is clear that the quartet phase is never energetically 
favorable to the macroscopically phase separated state. 
The fit to (|l4| ) is shown as the solid line in Fig. g. 

Because the four-body energy per particle is higher 
than that of the phase-separated state, we see no need 
to consider gases of larger clusters as potential ground 
states. In the limit of large particle number iV, the energy 
per electron of a finite bound droplet scales as 

£N^£N=oo+a/^/N + ---, (15) 

where £M=oa is the bulk energy per particle of the phase- 
separated state and a > is the surface correction. Since 
the surface correction is positive, the energies of large 
clusters approach the phase-separated energy monotoni- 
cally from above with increasing droplet size, and £n > 
£n=oo for sufficiently large N . Clearly this inequality 
is violated for A^ = 2, but we have shown it holds for 
iV = 4. Assuming the energy per particle of clusters is 
a reasonably smooth function of cluster size, the phase- 
separated state will have lower energy per electron than 
all droplets with four or more particles. Thus the zero- 
electron-density phase diagram consists only of the three 
phases described above. 

At low but non-zero electron densities, both the phase- 
separated and s-wave pair states remain robust. How- 
ever, the phase of free electrons at J < 2 is unstable to 
higher-angular-momentum pairings similar to the Kohn- 
Luttinger effect in three dimensions ||l3,nq|. Using ex- 
pressions (20) and (21) in Ref. |13], we find the phase dia- 
gram shown in Fig. 0, which differs substantially from the 



FIG. 3. Phase diagram for low electron density. The 
phase-separation boundary is estimated from the two-body 
solution on finite systems, and the other lines are calculated 
from expressions in Ref. ||l5|. All phase boundaries become 
less accurate with increasing electron density per site n. 



phase diagram in that Reference. As the electron density 
increases from zero, the phase-separated and s-wave pair- 
ing states persist, but the free-electron phase is unstable 
to p-wave pairing at small J and (1^2 _y2 pairing at larger 
J. The energy gaps in these higher-angular-momentum 
pairing states increase from zero with increasing density, 
but remain extremely small at the low densities where the 
effect is valid. For example, at electron density n = 0.1, 
all p and d-wave energy gaps are less than 10~^ in units 
of i. 

In summary, we have calculated the complete phase 
diagram of the two-dimensional t-J model in the limit of 
low electron density. The phase-separated state at large 
J evaporates into a gas of s-wave pairs as the interaction 
strength is reduced — Gases of larger finite clusters are 
never stabilized. The s-wave pairs unbind as J is reduced 
further, but the resulting Fermi liquid is unstable first to 
d-wave pairing and finally to p-wave pairing. 
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